Regular packings on periodic lattices 
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We investigate the problem of packing identical hard objects on regular lattices in d dimensions. 
Restricting configuration space to parallel alignment of the objects, we study the densest packing 
at a given aspect ratio X. For rectangles and ellipses on the square lattice as well as for biaxial 
ellipsoids on a simple cubic lattice, we calculate the maximum packing fraction ifid{X). It is proved 
to be continuous with an infinite number of singular points X™'", X™^'' u = Q, ±1, ±2, .... In 
two dimensions, all maxima have the same height, whereas there is a unique global maximum for 
the case of ellipsoids. The form of ifid{X) is discussed in the context of geometrical frustration 
effects, transitions in the contact numbers and number theoretical properties. Implications and 
generalizations for more general packing problems are outlined. 

PACS numbers: 61.50.Ah, 82.70.Dd, 64.70.kt 



The question of how densely objects can fill a vol- 
ume has attracted both, mathematicians and physicists, 
for centuries. One famous problem is that of pack- 
ing spheres. In 1611, Kepler conjectured that the face- 
centered cubic (fee) and hexagonal close packed (hep) lat- 
tice configurations of identical spheres yield the highest 
packing fraction (p'^^i^ — 7r/-\/l8 = 0.7404. Gaufi could 
show in 1831 that these are the optimal periodic packings 
of spheres, but only very recently it was proved that they 
are optimal within all possible arrangements [1 . Even for 
disks in the plane, the corresponding proof of optimality 
of the hexagonal packing with ipf^^ — ■k/\/V1 = 0.9069 
was only found in 1943 P]. Apart from its theoretical 
attraction along with its relation to coding theory [4], 
packing is a problem of practical relevance. Not only 
have practitioners long known that a densest packing 
of oranges or cannon balls can be achieved via hexag- 
onal layering but, more recently, packing problems have 
received substantial attention in engineering and opera- 
tions research as problems of optimizing yields in produc- 
tion or minimizing leakage currents in integrated circuits 
(see, e.g., Ref. [3]). 

In physics, periodic packings [IHS] are relevant for de- 
scribing and understanding crystalline materials. In con- 
trast, random close packings [7], i.e., maxima of the pack- 
ing fraction under some local dynamics starting from 
loosely packed configurations, have been used to model 
glasses [5] and granular materials [51. For spheres in 
3d, random close packing leads to a packing fraction 
'^'rcp ~ 0.64, significantly below 'p'^^- The hard ob- 
jects considered in such packings need not be spheres, 
but can be more general convex bodies. Although re- 
cently there has been extensive numerical work using 
techniques from dynamic programming and heuristic op- 
timization, complemented by experiments, for studying 
periodic packings [TO] or random close packing |11| for 
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non-spherical objects, there is a lack of analytical under- 
standing of these problems. For random close packing, 
it has been observed that the packing fraction increases 
over (/jJIqp ~ 0.64 as spheres are replaced by ellipsoids, 
and might even approach (ysj^^^ ^ 7r/-\/l8 = 0.7404 in 
some cases [11,. Concerning periodic packings, an affine 
transformation maps the fcc/hcp sphere packing to a pe- 
riodic lattice packing of identically aligned ellipsoids with 
maximum packing fraction (^jj^^ — '"'/'v/lS. That non- 
parallel arrangement of ellipsoids of revolution may ex- 
ceed ttZ-v/IS has been predicted in Refs. [HI [T3]. Such 
super-dense packings of ellipsoids were studied recently 
in more detail [T31 [H] . Particularly, it has been shown 
that LP = 0.7707 for aU aspect ratios X > [T4j. 

We make progress in the analytical understanding of 
the problem of packings of non-spherical bodies by taking 
a complementary approach. Instead of finding the lattice 
structure that maximizes the packing fraction for a given 
type K. of objects, we start out from a fixed Bravais lattice 
A and attach a body JC of the same shape and orientation 
u to each lattice site (at its center of mass, say). We then 
determine the maximum packing fraction as a function 
of IC, i.e., as a function of the parameters characterizing 
its shape and orientation. To the best of our knowl- 
edge, this problem has not been studied before. Our ap- 
proach may contribute to describing, for instance, plastic 
crystals, i.e., lattices with a molecule fixed at each site. 
In particular, aromatic molecules can be approximately 
described by hard ellipsoids. Similarly, applications are 
envisaged in operations research and manufacturing. Fi- 
nally, insight into the frustration effects generated by the 
competing length scales of K, and A could contribute to 
the understanding of packings without a pre-determined 
lattice structure. 

Consider a class of identical d-dimensional convex bod- 
ies K. whose shape depends merely on their "length" I 
and "width" w, and consequently are characterized by 
a single parameter X — l/w, the aspect ratio. As a 
general example one might think of a d-dimensional el- 
lipsoid of revolution. Fixing the aspect ratio X and ori- 
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FIG. 1. (Color online) Maximum packing configurations 
of rectangles with aspect ratios X — 4 (green, bottom), 
X = XJi'^" = 5 (blue, middle), and X = 7 (red, top), re- 



spectively. Lattice lines L 
are indicated with dashed lines. 
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entation u), proportional rescaling of the bodies allows to 
reach the maximum packing fraction without overlaps, 
ipd{X,ijj). This fraction varies with w, and we are in- 
terested in the maximum packing fraction irrespective 
of orientation, ifd{X) = max^^ <^)- The maximum 

(fd{X) is continuous as a function of X. Here, we only 
outline the idea of the rigorous proof [IB]. Let us as- 
sume that (pd{X) is discontinuous at some Xq = Iq/wq 
where, e.g. it jumps from ip^ to ip+ > tp^, with 
ip± = Vnnipd{XQ ± e). The convex bodies at ip_ and 

are characterized by (?_,?«_) and (Z+,w+), respec- 
tively. Both pairs differ from each other, as (y5_ =/= ip^. 
Of course, it is l+fw^ — l^/w^ — Xq. Now, starting 
from the configuration at we continuously decrease 
the length of the hard objects. Consequently, without 
change of orientation, both the aspect ratio X and the 
corresponding packing fraction ipd{X) decrease continu- 
ously from Xq and (p+ = ipd^Xo), respectively. Below 
but arbitrarily close to Xq, (pd{X) must be arbitrarily 
close to due to its continuity. On the other hand, 
it is PdiX) > (pd{X) for all X < Xq, since ipd{X) is the 
maximum packing fraction by definition. Therefore, even 
if (pd{X) — <PdiX) holds (instead of >) for all X < Xq, 
we get — lim ipd{XQ — e) — (^-|_. This contradicts the 

original assumption > ip^. Consequently, ipd must 
be continuous. 

We now turn to the calculation of Pd{X) for specific 
hard objects. As an example in two dimensions (2d), we 
study a square lattice with lattice constant a — 1. Con- 
sider first the case of rectangles of length I and width 
w. Imagine two identical rectangles with common direc- 
tion e = (cos sin w), of their long side, attached with 
their centers to lattice sites (0,0) and Hjk = {j,k), re- 
spectively. In the following, we assume that j > and 
A; > are coprime integers, i.e., they do not have a com- 
mon divisor other than 1. For fixed aspect ratio X, it is 
obvious that the rectangles will attain maximum volume 
V2{l,w) if they touch each other and line up precisely 
along their short or long sides, cL Fig. [l] Combining this 



FIG. 2. (Color online) Maximum packing fraction ip2{X) for 
rectangles (dashed line, right scale) and for ellipses (solid line, 
left scale). Inset: orientation angle cj as a function of X for 
rectangles (dashed line) and ellipses (solid line). The crosses 
correspond to the positions of the maxima of (p2{X) for el- 
lipses. 



and the periodicity of the packing, it is straightforward 
to prove that e must be parallel to Rjk [H]. In other 
words, for given X maximum packing fractions will al- 
ways occur for "rational" orientations ui = arctan(fc/j) of 
the rectangles. Maximal packings for specific X can be 
constructed using the concept of lattice lines. The line 
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■j-j, through the origin is defined by the lattice vector 
Rjk = {j,k)- The distance of adjacent lattice sites on 



'^p} equals Ijk = \/ p + k^- The square lattice can be 



decomposed into a set of parallel lattice lines , a — 0, 
±1, . . ., of distance Wj^, where Wjk Ijk = vq = 1 (cf. the 
dashed lines in Fig. [l|). Choosing I = Ijk and w = Wjk, 



+ k'^ leads to a perfect 
= Ijk Wjk = 1, for all co- 



i.e., X = Xfl^- = Ijk/wjk = j 
tihng with i^^ax = ((52(^jT'') = 
prime pairs (j, k), cf. the maxima at ip2 — 1 in the lower 
part of the main panel of Fig. [2j The pairs (j, k) can be 
ordered such that X^^f < X^^"" where X^^"" = Xf'^H , 
v — 0, 1, 2, . . .. Since e must be parallel to R^fc, the 
maximum packing for X < Xp,^^ and X > X^'^^ is ob- 
tained by decreasing / below Ijk keeping w = Wjk, and 
decreasing w below Wjk keeping I = Ijk, respectively (see 
Fig. [T]). Consequently, 
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follow from the 
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j.+ik^+i - ^- V>2{X) is shown in Fig. [2j 
together with the optimal orientation uj{X) in the inset. 
We now turn to the case of packing ellipses on the 
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square lattice. A naive approach would be to inscribe 
them into the rectangles considered above. The resulting 
packing fraction of ellipses is then just 7r/4 that of the 
rectangles. In reality, however, maximally packed ellipses 
do not, in general, touch each other "head" to "tail", nor 
are they oriented parallel to the lattice lines, cf. Fig. [3] 
In contrast to the highly degenerate case of packing rect- 
angles which touch along whole line segments, packings 
of general, smooth convex bodies are characterized by K 
contact points per body of which, due to inversion sym- 
metry, only K/2 are independent. The three parameters 
describing an ellipse (two half axes and the orientation 
angle) are under-determined in the generic case of if = 4 
contact points (resulting in K/2 = 2 equations), yielding 
a continuum of solutions as a function of X. Non-generic 
is the case oi K — Q contacts, leading to a discrete set 
of maxima in ip2{X). For this situation, put one ellipse 
at the origin, such that the sites of the other ellipses are 
at ±(^i, TTii), i = 1, 2 and ±(j, k) — ±(Zi -I- h^mi + 7712), 
cf. Fig. [3| Note that < < j, < m,; < fc. Then, 
the three corresponding contact vectors Ci — ^{li,mi), 
i = 1, 2 and C3 = ^{j, k) = Ci + uniquely determine 
the three coefficients a, 6, and c in the ellipse equation 
ax'^+2bxy+cy'^ = 1. This allows to determine the lengths 
of the half axes and thus the aspect ratio to be 

Xffe- = {a+ + .^al+al)/Vi, (2) 

tti = Ci • £jCi + C2 • £iC2 + Ci • eiC2, 

where Sq — and £± = ^ . The correspond- 

ing packing fraction is ipn^^x = 'P2{Xjl^^) — 7r/\/T2 for 
all {j,k), identical to the packing fraction of hep disks. 
In fact, each such maximal ellipse packing can be con- 
tinuously deformed via an affine transformation into a 
packing of disks. On increasing (decreasing) X from 
Xp.^^, the shortest (longest) contact vector disappears, 
and the remaining four contacts allow to determine the 
coefficients a, b and c, and therefore ip2{X) and w(X) 
in between the maxima X^^^ as a function of X in a 
closed-form expression. 

The result for ip2{X) and uj{X) is displayed in Fig. [2] 
Since (p2{l/X) = (p2{X), only the regime X > 1 is shown. 
The maximum packing fraction is singular at X^^^ for 
all ly. The orientation uj{X) is discontinuous at those 
X™™ at which (p'2{X) is discontinuous and at those X^l^^ 
which are degenerate, such as X^'^^ = X^'^^ (cf. Fig. |2]) . 
At these points, there are two degenerate maximal pack- 
ings with the same packing fraction and aspect ratio, 
but different orientations. The global maximum value of 
fmax = 1 and (/9,nax = TT / y/l2 for rcctaugles and ellipses, 
respectively, is attained for an infinite number of pack- 
ings, uniquely labeled by {j,k). From Fig. |2]it appears 
plausible that lim ip2iX) — (pmax, which indeed can be 

proved [IB] . 

The relation of the contact points can be understood 
from a number-theoretical point of view. Note that the 



FIG. 3. (Color online) Maximum packing configurations of 
ellipses with aspect ratios X — 3 (green, bottom), XJI^" = 

\/29-f 8^13/^3 = 4.4 (blue, middle), and X = 6 (red, top). 
Lattice lines for {j, k) = (2, 1) (dashed lines) and contact 
points (crosses). 

centers {li,mi), i = 1, 2 of two ellipses touching the cen- 
tral one at (0, 0) also define lattice lines with direc- 
tion {li,mi). These are the directions closest to that of 

-^jfe' provided that U and irii are coprime and < < j, 
< rui < k. In mathematical terms, this means that 
nii/li, z = 1, 2, are given by the best principal and best 
intermediate rational approximant [17] of k/j. They fol- 
low from the finite continued fraction expansion of k/j, 

k/j = ao + + l/[a2 + • • ■ + l/K-i + l/a™] ••■]], (3) 

where a^, i — l,...,n (a„ > 2) are positive inte- 
gers that are uniquely determined by k/j. Then, it is 
h = Sri-ii ™i = where the best principal approx- 

imant r„_i/s„_i follows from Eq. ^ for a„ = 00, and 
I2 = Sn,a„-i^ ™2 = fn,an-i foUows analogously from 
the best intermediate approximant rn.a„-i/ Sn^a^-i ob- 
tained from Eq. ([3| replacing a„ by a„ — 1. Since coding 
problems are strongly linked to number theory [3], these 
results also promise insight into the connection between 
packing and coding problems. 

Finally, we have investigated ellipsoids of revolution 
on a simple cubic lattice. Analytically, it is possible to 
proceed in a similar fashion as for the ellipses. The re- 
sulting eighth-order polynomial in can only be solved 
numerically, however, and the intermediate expressions 
are rather unwieldy. Therefore, we instead determined 
ipz{X) numerically by a downhill-simplex minimization 
algorithm, the result of which is shown in Fig. [4] it agrees 
with that determined earlier in Ref. [TH] and, as expected, 
shows continuity, too. Similar to the results in 2d, the 
derivative ip'^{X) seems to be discontinuous at a series 
of maxima at X™^'^. It appears to be discontinuous at 
some, but not all, minima X™™. The symmetry between 
1/X and X valid in 2d, however, is lost. Most strik- 
ingly, the global maximum (yJmax = 7r/-\/l8 for ellipsoids 
appears to be attained only for the single packing frac- 
tion X'^i^ — 1/2, whereas in 2d there was a countable 
infinity of degenerate maxima. This maximum corre- 
sponds to the highly non-generic case of each ellipsoid 
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FIG. 4. (Color online) Maximum packing fraction for biaxial 
ellipsoids on an sc lattice on a logarithmic X scale. 

touching 12 neighbors. Consequently, an affine trans- 
formation can be apphed to map this pattern to closest 
packing of hard spheres resulting in an fee or hep lat- 
tice. A second prominent maximum occurs at X^'^'^ = 2, 
where 8 contact points occur. The corresponding trans- 
formed hard sphere packing yields a bcc lattice. It can 
be shown that (fisiX) > (/Smin = 7r/6 = (p^^X = 1), i.e., 
the packing fraction of hard spheres on an sc lattice is 
a lower bound for tp^(X) |16| . From the numerical re- 
sults in Fig. |4] we conjecture that, analogous to the 2d 
case, (psiX) -> <Pmax — 7r/\/l8 for X oo and X — > 0, 
respectively. 

To conclude, the maximum packing fraction ipd{X) of 
parallely aligned convex objects /C characterized by a sin- 
gle aspect ratio exhibits universal features that appear 
to be independent of /C, the underlying Bravais lattice A, 
and even its dimension d. In particular, ipdiX) has been 
very generally proved to be continuous. In fact, this proof 
can even be extended to the case of convex bodies char- 
acterized by an arbitrary number of aspect ratios. Fur- 
thermore, as shown for the case of rectangles and ellipses 
on the square lattice as well as for biaxial ellipsoids on 
the sc lattice, there is an infinite number of local maxima 
and minima at which Lpd{X) is singular. The singular- 
ities at the minima and at certain, degenerate maxima 
(see Fig. [2] for the case of ellipses) are correlated to the 
discontinuities in the orientation of /C. For the studied 
cases, we find that ipd{X) converges to its global maxi- 
mum for A — oo as well as for A' ^ 0. While we were 
only able to prove this rigorously for the case of rectan- 



gles and ellipses, we believe that this property holds far 
more generally, implying that convex hard objects, on 
average, pack much better if they become more oblate or 
prolate. 

On the other hand, there are also significant differences 
between the systems studied in two and three dimensions. 
For rectangles and ellipses, the global maximum pack- 
ing fraction y^max = 1 (rectangles) and (praax = 7r/-\/l2 
(ellipses) is attained for an infinite number of discrete 
aspect ratios XJ^^^, uniquely labeled by pairs (j,fc) of 
coprime integers. On the contrary, for symmetric ellip- 
soids with 1/100 < A < 100, <f3{X) takes its maximal 
height (^niax — 7r/\/l8 at the single value A^f^ = 1/2, 
only. This qualitative difference can be understood as fol- 
lows. Consider, for instance, a d-dimensional symmetric 
ellipsoid which depends on d + 1 parameters. In a pack- 
ing, K contacts lead to K/2 equations. In the generic 
case of K/2 = d, the system is under-determined and (pd 
can be found as a function of the aspect ratio A. For the 
non-generic case K/2 = d -I- 1, there is always a solution 
corresponding to the local maxima of ipd{X) at A™'^'^. It 
appears likely that the competing point symmetries of /C 
and A are responsible for the non-equal heights of these 
maxima for c? = 3 and K = 8. It is conceivable that this 
extra frustration might be relieved by considering convex 
hard objects characterized by three length scales, possibly 
leading again to an infinity of equal-height maxima. The 
even more non-generic situation K/2 > d+l as realized, 
e.g., in the global maximum (/^max = 7r/\/r8 for our 3d el- 
lipsoids with K — 12, corresponds to an over-determined 
set of equations such that, at most, only very few solu- 
tions can be expected. It is worthwhile to point out that 
our results for packing on fixed lattices should be closely 
related to "continuum packing with a fixed number of 
contacts since the latter involves geometric frustration 
as well. 

Of course, it might be a challenge to study packings 
of the considered type on different lattices. Even richer 
behavior is expected on weakening the condition of par- 
allel alignment, paving the way for the occurrence of su- 
perdense packings in analogy to those recently found for 
ellipsoids in the 3d continuum [TH [T5] . 
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